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Abstract 

The problem of statistical description of classical matter, represented by a N-body 
system of relativistic point particles falling into a black hole, is investigated, adopt- 
ing a classical framework. A covariant microscopic statistical description is devel- 
oped for the classical particle system by introducing the definition of observer- 
related N-particle microscopic distribution function and by suitably extending to 
general relativity the concept of Gibbs microscopic entropy. The corresponding en- 
tropy production rate is calculated taking into account the presence of an event 
horizon and proven to obey, for N » 1, an exact microscopic H-theorem which 
holds for arbitrary geometries. 
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1 Introduction 

Since its inception [H Q, 0, 0, II] the thermodynamical interpretation of black 
holes (BH) has been the subject of debate f^. Indeed, the mathematical anal- 
ogy between the laws of thermodynamics and black hole physics following 
from classical general relativity still escapes a complete and satisfactory in- 
terpretation. In particular it is not yet clear whether this analogy is merely 
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formal or leads to an actual identification of physical quantities belonging to 
apparently unrelated frameworks. The difficulty is related to the very concept 
of entropy usually adopted in BH theory, based on Boltzmann entropy, which 
is determined by the number W of microscopic complexions compatible with 
the macroscopic state of a physical system 



Sbh = K\nW. 



Indeed S^h does not rely on a true statistical description of physical systems, 
but only on the classification of the microstates of the system, quantal or 
classical. As a consequence, the evaluation of Sbh requires the knowledge of 
the internal structure of the BH, a result which obviously cannot be achieved 
in the context of a purely classical description of BH. Therefore the evaluation 
of Sbh requires a consistent formulation of quantum theory of gravitation and 
matter This can be based, for example, on string theory_0L and can be 

conveniently tested in the framework of semiclassical gravity js], |9| . 

A basic difficulty of quantum theories founded on the concept of Boltzmann 
entropy, is that up not now they have not leaved up to their expectations since 
they have not yet achieved their primary goal, i.e., the rigorous proof of an 
H-theorem and the full consistency with the second law of thermodynamics 
6Sbh > 0. Indeed, estimates of the Boltzmann entropy based quantum theory 
of gravitation 0, |^ and yielding Sbh = being A the area of the event 

horizon, are inconsistent from this viewpoint, since as a consequence of the 
BH radiation effect 0] the radius of the BH may actually decrease. Hence, the 
Boltzmann entropy Sbh cannot be interpreted, in a proper sense, as a physical 
entropy of the BH. To resolve this difficulty a modified constitutive equation 
for the entropy was postulated [l[ 0], in order to include the contribution of 
the matter in the BH exterior, by setting 



S' — S -\- Si 



bh, 



{S' denoting the so-called Bekenstein entropy) where 5* represents the cor- 
rection carried by the matter outside the BH (notice, however, that also S 
cannot be interpreted as entropy). As a consequence a generalized second law 
5S' > was proposed [H, 01 which can be viewed as nothing more than the 
ordinary second law of thermodynamics applied to a system containing a BH. 
However, the precise definition and underlying statistical basis both for S and 
potentially also of Sbh remain obscure. Thus a fundamental problem appears 
their precise estimates based on suitable microscopic models. 

On the other hand, if one regards the BH as a classical object in the space- 
time continuum, provided the surrounding falling matter can be assumed as 
formed by a suitably large number of particles, the estimate of the BH entropy 
should be achievable directly in the context of classical statistical mechanics, 
by adopting the customary concept of statistical entropy, i.e., Gibbs entropy. 
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In contrast to the Boltzmann entropy, this is based on a statistical description 
of physical systems and is defined in terms of the probability distribution of 
the observable microstates of the system. In fact, as is well-known, its defini- 
tion coincides with the axiomatic definition of Shannon entropy, yielding the 
measure of ignorance on statistical ensembles. A first result of this type been 
presented in previous work by Nicolini and Tessarotto (hereon denoted as 
Ref.I), where a covariant kinetic theory in the presence of an event horizon was 
developed for classical matter in the BH exterior, treated as a N-body system 
(with iV 1) of classical point particles undergoing, before capture, a purely 
Hamiltonian dynamics. By introducing a suitable definition for the relativistic 
kinetic entropy of infalling classical matter described as an ensemble of point 
particles, an H-Theorem was reached, tanking into account the presence of the 
black hole event horizon. 

The goal of the present paper is to extend the results of Ref.I by developing 
a microscopic (N-body) statistical description of classical matter in the pres- 
ence of an event horizon. For this purpose we intend to evaluate the Gibbs 
(microscopic) entropy of classical matter falling into the BH event horizon. In 
particular, we intend to prove the validity of an exact H-theorem holding for 
Gibbs entropy, provided the infalling matter can be described as a suitably 
large classical system (A^ ^ 1) of point (neutral or charged) particles forming 
plasma or a rarefied gas which interact mutually only via a mean Hamiltonian 
field. 

The scheme of the presentation is as follows. In Sec. II, first the basic as- 
sumptions of the theory are introduced, which include the concept of particle 
capture domain, assumed to be defined by a smooth hypersurface located 
suitably close to the event horizon {capture surface). Second, a covariant mi- 
croscopic statistical description is formulated for a N-body system of classical 
point particles referred to an arbitrary observer. In particular, in the exterior 
domain of the BH (outside the capture surface) the particle system is assumed 
to obey a Hamiltonian dynamics. As a consequence, in the same domain the 
Liouville theorem for the corresponding phase-fiow is readily recovered and, 
based on the concept of observer-related N-particle microscopic distribution 
function and probability density, the relativistic N-body Liouville equation 
is determined. Finally suitable boundary conditions are introduced for the 
microscopic distribution function on the capture surface. 

In Sec. HI the definitions of observer-related relativistic Gibbs and kinetic 
entropies for the infalling classical matter are introduced and the correspond- 
ing entropy production rates are determined. The relationship between the 
two entropies is proven to hold in an elementary way for classical systems 
of weakly interacting point particles. As a consequence the obvious physical 
interpretation of the contributions of the entropy production rate are pointed 
out. 
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Finally, in Sec. IV, an H-theorem is recovered for the Gibbs entropy, extending 
the result of Ref.I. The result is proven to hold for BH characterized by an 
event horizon of arbitrary shape and size. In particular, the result is proven 
to apply, in contrast to Boltzmann entropy in quantum theory of gravitation 
j^, also to the case of matter falling in a BH with contracting event horizons, 
such as BH implosions or slow contractions. 



2 The N-body covariant microscopic statistical description 



In this section we introduce the basic framework of the theory. The assump- 
tions we are going to make deal with the treatment of non-isolated black hole 
formed by the collapse of a star and surrounded by matter. The BH, together 
with its associated event horizon, and of the matter surrounding the BH and 
falling toward the event horizon are all assumed as classical. It is sufficient for 
our purposes to assume the following hypotheses: 

I) Infalling particles capture is due to the redshift phenomenon occurring near 
the event horizon for an arbitrary observer located far from the BH (for ex- 
ample in a region where space time is asymptotically fiat fl^I ). As a result 
particles sufficiently close to the event horizon effectively disappear to the ob- 
server (Assumption 1). In the sequel we shall assume that all particle capture 
events occurs in a subdomain, to be identified with a surface 7 of the space- 
time [capture surface), localized infinitesimally close to the event horizon. 

II) The total energy of infalling particles is finite in such a way that local 
distortions of space-time are negligible (Assumption 2). 

The matter outside the BH is described by a system of ^ 1 identical 
classical particles to be referred to an arbitrary observer O. If we assume that 
the system is Hamiltonian, its evolution is well known and results uniquely 
determined by the classical equations of motion, defined with respect to the 
observer O. To this purpose let us choose O, without loss of generality, in a 
region where space time is (asymptotically) fiat, endowing it with the proper 
time r, where r is assumed to span the set / C (observer's time axis). 
Without loss of generality we can assume that the particles are points, i.e., they 
are described by the 1-particle canonical states Xj (with i = 1, N) spanning the 
8— dimensional phase space Fj, where Xj = ir^.pi^). The analogous treatment 
of particles having higher degree of freedom is straightforward. Therefore, the 
evolution of the system, described in terms of the N-body canonical state x = 
{xi, ...,X7v} , is determined by a suitable relativistic Hamiltonian H = H (x), 
where each canonical 1-particle state Xj (z = 1, N) results parameterized in 
terms of the i—ih particle world line arc length Si (see ll2|]). As a consequence, 
requiring that = Sj(r) results a strictly monotonic function it follows that. 



4 



the particle state can be also parameterized in terms of the observer's time r. 
Therefore, the particle states are determined by the canonical equations: 



dsj^i) dxi{si{T)) _ ds(i) 
dr ds(i) dr 



H 



where 



/(x(j)), H is the canonical Poisson equation 



(2) 



g/(x(.)) dH _ a/(x(,)) dH 
drt dp(i)^, dp. 



(3) 



and the previous initial- value problem admits by assumption a unique solution 
defining a C^'^\r^ x 7)-solution. It follow that the phase- flow defined by the 

mapping Xo = {xoi, Xoat} ^ x(r) = {xi(si(r)), , XAr(sAr(r))} satisfies 

a Liouville theorem, namely ^ 
canonical measure defining the p 



ax(T) 



(see Appendix) and therefore the 
n dyii is conserved, 

i=l,N 



(9x0 

lase space volume dx 



since dx = dx„ 



Next, to obtain a microscopic statistical description for the N-body system we 
introduce the N—body microscopic distribution function for the observer O, 

pI^\x) = ^^(x) n - n H^/v^ - 1) = iv-pg"(x) (4) 



i=l,N 



-AM 



defined on the extended N— particle phase-space — f] T^'* (where 

i=l,N 

r^'\ for i = 1,A^, is the 8-dimensional phase space of the i-th particle).. 
Here p^^^ (x) denotes the conventional microscopic distribution function in 
the 8A^— dimensional phase space, while Pq^^ (x) is the corresponding A^-body 
microscopic probability density satisfying the normalization 



/ 



cixpjf^(x) 



(5) 



Notice that the Dirac deltas introduced above must be intended as physical 
realizability equations. In particular the condition placed on the arc lengths Sj 
implies that the i-th particle of the system is parameterized with respect to 
Sj(r), i.e., it results functionally dependent on the proper time of the observer; 
instead the constraints placed on 4- velocity implies that -Uj^ must belong to 
the hypersurface di of equation ^Ju^jV^ = 1, and hence is a tangent vector 
to a timehke geodesic. In the sequel we adopt also the notation 



p<^)(x) 



p(^)(x) 



(6) 
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to denote p(^^(x) evaluated on the intersection of the hypersurfaces Si, ..... S^. 
The event horizon of a classical BH is defined by the hypersurface r^f, identical 
for all particles and specified by the equation 

Rm = th, (7) 

where for m = 1,N, Rm denotes the a suitable curvilinear coordinate for 
the m-th particle and coincides with the radial coordinate in the spherically 
symmetric case. According to a classical point of view, let us now assume that 
the particle capture surface be defined by the surface 7 of equation 

Rm{Sm) ^r^{Sm), (8) 

where = (1 + e)rjy, while e > is an infinitesimal which may depend on ex- 
plicit or hidden parameters (for example, e might depend on the detector used 
by the observer). The presence of the BH event horizon is taken into account 
by defining suitable boundary conditions for the kinetic distribution function 
on the hypersurface 7. For this purpose we distinguish between incoming and 

outgoing distributions on 7 with respect to the i-ih. particle, p[^'"-^'*'* (x) and 

7 

corresponding respectively to > and UaUf < 0, where Ua 

7 

is a locally radial outward 4— vector. Therefore the boundary conditions on 7 
are specified as follows 



pg^^(-'^)(x)i ^pW(x) n ks,-s,{t)) n 5(^^-1), (9) 

^ i=l,Af i=l,Af 



= 0. 
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As previously anticipated, these boundary conditions do not actually require 
the detailed physical model for the particle loss mechanism, since all particles 
arc assumed to be captured on the same hyper-surface 7, independent of their 
state. This provides a classical loss model for the BH event horizon. 



3 Covariant Liouville equation 



It is now immediate to determine the covariant Liouville equation which in 
the external domains, i.e. outside the event horizon, advances in time the 
microscopic distribution function Pq\^ with respect to the observer O. 
Thanks to Liouville theorem (i.e., phase-space volume conservation in in 
the sense indicated above) and invoking as usual the axiom of conservation 
of probabihty for classical systems, it follows that p^^\x.) must satisfies the 



6 



differential Liouville equation 



ds^(d4dp^2M + ^^EM]=o (10) 
dr \ dsi c?rf dsi dpi^ J 

where by assumption 

^>0 (11) 
dr 

is made for all z = 1, A^, and the summation is understood over repeated 
indexes. This equation resumes the conservation of the probability in the rela- 
tivistic phase space in the domain external to the event horizon. Invoking the 
Hamiltonian dynamics (2), the kinetic equation takes the conservative form 



dsi 
dr 



p(^)(x),if 



0. (12) 



Let us now introduce the assumption, holding for collisionless of weakly inter- 
acting particles, either neutral or charged, that the Hamiltonian H (x) can be 
expressed in the form 

H{^)= Y: H„ (13) 

i=l,N 

where Hi = Hi (xj) . This implies that the particles of the A^-body system 
interact mutually only via a mean-field Hamiltonian force. It follows that the 
Liouville equation (10) admits the factorized solution 

p(^)(x) = iV n P(xO (14) 

i=l,N 

where p(xj) = p^^\xi)/N is the kinetic (1-particle) probability density and 
p(xj) is the related kinetic distribution function, which manifestly obeys the 
covariant kinetic equation of the form 0, |2^ 



dsii 



dr 



p«(x(i)),iJ(J ^=0. (15) 



Notice that this equation is independent of N, the number of particles, to be 
assumed in the sequel as finite. 



4 Observer-related relativistic Gibbs and kinetic entropy 



Let us now introduce the definition for the microscopic entropy S'(p''^^) appro- 
priate in the context of the present covariant theory. The definition follows by 
straightforward generalization of the relativistic kinetic entropy llOU a nd the 
concept of Gibbs entropy in non-relativistic [3, 3] and relativistic |l6l. 17, 3| 
statistical mechanics. Thus the concept of Gibbs microscopic entropy can be 



7 



defined in analogy to Ref. I, with respect to an observer endowed with proper 
time T 



/ d^{s) n s{sj-sj{r)) n (le) 

pi, j=l,N k=l,N 

Here the notation is as follows: dx(s) = H d'^ji^j): where for each particle 

j=l,JV 

j — 1,N the state vectors Xj are parameterized with respect to Sj, with Si 
denoting the s-particle arc length. Finally, P denotes the principal value of 
the integral introduced in order to exclude from the integration domain the 
subset in which the distribution function vanishes. It is immediate to obtain 
the relationship between Gibbs and kinetic entropy, S{p^^^), previously defined 
in Ref. I and given by the equation 



Sip^""^) = - / c^xi(s)5(si -si(r))5(y^- iyi)lnp«, (17) 
r 

where F denotes the 1-particle phase space. In fact the condition of factoriza- 
tion (14) implies immediately 

5(p(^)) = Ar5(pW). (18) 

The kinetic entropy can also be written in the equivalent way 

5(p(i)) ^-Pj dMs)S{s^ - Si(t))pS'Hx) lnp(^)(x), (19) 
r 

where p^\'x(s)) reads 

p?{^{s)) = e{n{s^) - r,{s^))5{^/^,-l)p^'\^{s)), (20) 
where © denotes the strong Heaviside function 

1 for a> 

e(a) = <^ (21) 
for a< 0. 

Equations (18) and (19) allow to determine immediately the entropy produc- 
tion rate associated to the Gibbs entropy, which reads 



ci5(p(^)) 



dr 



^-N I dh,d'p^FrrJ{r,-r,,)p^'^]np^'^ = S, + S2, (22) 



where F is the subdomain of phase space in which riaUf < and F^r^ is the 
characteristic integrating factor 

_ dsiir) / dfi dr^i\ 
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It follows that ^^^^^and can be interpreted as the average of the entropy 

flux across the capture surface 7. Moreover, and denote respectively the 
contributions to entropy production rate 

S, = -NP J ^'rid^pi^^5 (n - lnp«, (24) 

r- 

S.^NPJ d'r.d'pi^^^S (n - p« lnp«. (25) 
r- 

Wc stress that here by construction ^^^^-^ > [see Eq.(ll)] and ^ < 0, ^ 
denoting the "radial" velocity of infalling matter on the surface 7, while there 
results < 0, = or > 0, being ^-the local velocity of the surface 7, re- 
spectively for contracting, stationary and expanding event horizons. However, 

the signs of and 5*2 are generally not defined, unless further assumptions 
are taken into account. Notice that, in analogy to the Bekenstein position (1), 

(Si and 5*2 denote the contributions to the entropy flux carried by incoming 
matter and by the BH due to the motion of the event horizon, therefore they 
can be at least qualitatively related in the following way: 



Si^S, (26) 

being S and Sbh the contributions to Bekenstein entropy (1). In particular, it 

interesting to remark that both ^land S2 result by construction proportional 
to A, the area of the event horizon, a conclusion which appears in qualitative 
agreement with estimate for the BH Boltzmann entropy given above for Sbh- 



5 H-theorem for the Gibbs entropy 



Let us now introduce the assumption that the total number of particles is 
finite, but suitable large {N ^ 1). In such a case, it is possible to determine 

the signs of Si, which contributes to the entropy production rate (22). Indeed 

it is possible to prove that there results 5*1 > while S'2 has not a definite sign. 
In addition, thanks to the results given in the previous section, in particular 
the relationship between the Gibbs and kinetic entropies, S{p'^^^) and S'(p(^^), 
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specified by Eq.(18), the following H-theorem holds for S{p^^^), 



^^.5. + S,>0. (27) 
ar 

Moreover, we notice that the support of the kinetic distribution function, i.e., 
the subset of in which the kinetic distribution function is non negative, 
results always compact. This condition direct consequence of the As- 

sumption 2 here considered, implying that the energy of the falling particles 
reaching the surface 7 cannot become infinite. 

Therefore, denoting Q the subset of F" in which the kinetic distribution func- 
tion p^^^ is non-zero, in the complementary set F" \ fi, the kinetic distribu- 
tion function p^^^ = Np^^^ (being p^^^ the kinetic probability density) results 
identically zero. Thanks to Assumption 2, it follows that such a set results 
necessarily bounded. Therefore, so that the following majorization holds 



1)1 



>pj dhd^p \Frr, \S{r- re) [Arp(i) - l] . (28) 
Thus, letting 

Ms = jdhd^p\Frr,\5{r -Te) (29) 
and imposing that A^ ^ 1 be sufficiently large to satisfy the inequality 

d5(p(i); 



dT 



> ^ = AT inf 1^ dhnoV;^^ | - M5 > (30) 



the thesis of the H-theorem is reached provided we assume inf 
0, a condition consistent with the requirement of a non isolated BH surrounded 
by matter. In the previous equation we have introduced the additional nota- 
tion 

P f dh,d'p,^^^S (n - ra) p(') = N [ dhnoV;", (31) 
J or dsi JQ 

r- 

being uq the number density. We stress that this result generalized the H- 
theorem given in Ref.I, since it applies also to Gibbs entropy. The result holds 
for classical BH having, in principle, arbitrary shape of the event horizon and 
even in the presence of a contracting event horizon (which might by produced, 
for example, by star implosions). The present theory appear therefore poten- 
tially relevant for a realistic detailed analysis of the BH thermodynamical 
properties. 
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6 Conclusions 



In this paper a macroscopic statistical description has been adopted for clas- 
sical matter around black holes. Matter in the immediate vicinities of a BH 
event horizon has been modelled by a weakly interacting relativistic gas Sn- 
Its dynamics results described by the relativistic Lioville equation, while the 
presence of the BH event horizon is taken into account by treating it as a 
classical absorbing porous wall. 

By assuming that Hamiltonian dynamics takes into account only mean field 
interactions between particles, the connection with the kinetic treatment of 
Ref.I can be immediately estabhshed. As a consequence, an H-theorem valid 
for the Liouville equation can be established on rigorous grounds which applies 
to every space time geometry and to the case of contracting horizon. 
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8 Appendix: relativistic N-body Liouville theorem 



Let us assume that the point particles constitute an isolated N-body system 
{assumption a) obeying the relativistic Hamiltonian equations of motion 

= X„ (32) 

asi 

Xi(sJ = Xoi, (33) 
which imphes that all the vector fields Xj (i = 1, N) are conservative, i.e., for 

■ X, = 0. (34) 

<9x(i) 

Introducing the parametrization in terms of the observer's time r and requiring 
that the functions Si = Si(r) are strictly monotonic, the equations of motion 
can be written in the symbolic form it follows 



dr dr dsj dr 



X,. (35) 
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Here x = {xi, XAr}^ = {yi, ,y8N} , X = {Xi, X^vj^ = {Yi, , Ysn} , 

where by assumption a the vector field X depends only on the (local or re- 
tarded) states of the particles forming the A^-body system. As a consequence, 
let us denote 



x(t) = x(s(t)) = {xi(si(t)), xjv(s;v(t))} 



(36) 



the solution of the initial value problem (32)-(33) and x(ro) = x(s(ro)) = 
{xi(si(to)), Xjv(siv(To))} = Xq the initial condition, where x(t) and x(to) 
denote the A^— body system states as seen by the observer O, respectively at 
times T, To- The previous assumptions for the phase- mapping 



Xo x(r) = x(x„. To, t) 
imply the following theorem: 



(37) 



8.0.1 THM. - Relativistic N-body Liouville theorem 
For arbitrary G F and Tq, r e / C M there results 

9x(t) 



1. 



(38) 



In fact, for A" > 1 the time derivative of the Jacobian 



9x(r) 
9xo 



reads 



d 

dr 



9x(r) 



<9xo 



E 



i=l,8Ar 



Hence by the chain rule 



= S I? E 

i=l,8Ar r=l,8Ar 



dsi 

d^ 

d_ 
d^ 



dxir) 



dr „j^^:,r dyrir) 



)) 



9X0 

9(y,(r),....,fi,...,y8iv(r)) 



9x(r) 



9xo 



^ dsj dYj 
i=MJV dr dyi 



9x(t) 



9xo 



9xo 



i=l,N 



dr dxi 



and thanks to the condition of conservation (34) 



(39) 



(40) 



d 

d^ 



ax(r) 



9xo 



0, 



which imphes the thesis, c.v.d. 



(41) 
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